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Introduction 

The  field  emission  current  fluctuation  method  for  determining  surface 
diffusion  coefficients  of  adsorbates  on  metal  surfaces,  proposed  by  the 
author  in  1973  [l,2]  has  been  applied  successfully  to  metals  which  can  be 
fabricated  into  clean  field  emitters.  It  has  recently  been  extended  to 
measurements  of  diffusion  anisotropy,  i.e.  a  determination  of  the  2-dimen¬ 
sional  diffusion  tensors  of  adsorbates  [3,4]  and  also  provides  the  compres¬ 
sibility  of  the  adlayer  [5]  and  information  on  correlation  lengths  and  their 
anisotropy  [3,4].  A  number  of  quite  unexpected  and  interesting  results 
have  been  obtained  so  far,  which  suggest  that  it  would  be  worthwhile  to 
extend  the  method  to  substrates  which  cannot  easily  be  fashioned  into  clean 
field  emitters.  The  recent  advent  of  scanning  tunneling  microscopy  (STM) 

[6]  suggests  some  extensions  of  the  method,  which  are  outlined  in  this  paper 
together  with  calculations  aimed  at  quantifying  them.  Finally  some  discus¬ 
sions  of  experimental  constraints  and  feasibility  criteria  are  given. 


Single  Site  Correlation  Function 

The  density-density  time  autocorrelation  function  evaluated  at  a  single 

site  provides  essentially  the  same  information  as  the  "patch  correlation 

function",  used  experimentally  to  date  [1,2].  Monte  Carlo  simulations  by 

Tringides  and  Gomer  [7]  use  single  site  correlation  functions  extensively 

to  obtain  D.  The  function  can  be  obtained  in  terms  of  t/x  where  the  relax- 

o 

ation  time  xQ  is  here  given  by 

xo  =  a2/ 4D  (1) 

a  is  a  dimension  corresponding  approximately  to  a  substrate  atomic  radius. 

If  the  diffusion  coefficient  D  is  expressed  as 
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D 


a  ve 


-E/kT 
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(2) 
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it  is  seen  from  Eq.  1  that  tq  takes  the  particularly  suggestive  form 

,  -E/kT, -1 

tq  =  (v  e  )  (3) 

where  v  is  an  effective  attempt  frequency  and  E  the  activation  energy. 

The  smallness  of  tq  for  this  case  brings  advantages  and  disadvantages. 

The  main  disadvantage  is  that  a  time  xQ  should  contain  100-200  individual 

correlation  measurements  or  time  boxes  At  so  that  At  ^  0.01  t  .  Ibis  requires 

either  that  only  small  values  of  D  be  investigated  or  that  rather  fast  electro- 

-92-1  2  -15  2 

nics  be  employed.  For  instance  if  D  =  10  cm  sec  and  a  =  10  cm 

_7 

=  2.5  x  10  sec.  At  'v.  2.5  nanosec.  While  this  may  be  possible  it  would 

-11  2  -1 

probably  be  necessary  to  work  at  temperatures  so  low  that  D  <  10  cm  sec 

in  order  to  increase  tq  and  At  accordingly.  The  advantage  of  small  tq  is  that 
it  may  minimize  drift  and  vibration  problems,  in  that  a  single  sequence  of 
200-400  channels  containing  one  correlation  measurement  can  be  carried  out 
over  a  given  site  without  interference.  Since  the  averaging  required  for  a 
meaningful  correlation  function  need  not  be  restricted  to  the  same  site  but 
can  be  carried  out  by  summing  measurements  over  equivalent  sites  slow  lateral 
drift  should  be  no  problem.  Vibrations  which  vary  the  tip-substrate  spacing 


can  also  be  removed  in  the  single  site  mode,  since  it  should  be  possible  to 


clip  the  signal  in  order  to  make  it  uniform.  All  that  is  required  is  that 
the  "occupied"  signal  differ  sufficiently  from  the  "empty"  signal.  In  any 
case  vibrational  periods  will  generally  be  «  x  . 

The  mean  square  fluctuation  is  obtained  as  usual  (2)  as  the  zero  time 
correlation  function.  To  convert  the  observed  value  to  a  number  fluctua¬ 
tion  function  it  is  only  necessary  to  divide  by  the  square  of  the  amplitude 
corresponding  to  the  "occupied"  signal.  It  is  perhaps  worth  pointing  out 
that  even  at  very  low  coverages  where  an  ad- atom  wanders  into  the  field  of  view 
only  occasionally  the  method  gives  the  correct  result,  since  the  time  averaging 


involved  in  obtaining  a  correlation  function  automatically  reduces  its 
amplitude  to  the  value  corresponding  to  the  particular  coverage  used. 

Ideally  the  latter  should  be  determined  independently.  It  may  happen, 
however,  that  adsorbate  coverage  is  very  low,  in  which  case  at  least  an 
estimate  may  be  obtained  from  f(o)  since  for  very  low  coverages  this 
quantity  is  equal  to  the  mean  particle  number  in  the  region  of  observation, 
i.e.  directly  gives  the  adatom  to  substrate  atom  ratio  [5]. 

Patch  and  Line  Correlation  Functions 

Although  the  single  site  method  offers  some  attractive  features  it  is 
limited  to  fairly  low  values  of  D  as  we  have  seen.  Thus  it  is  worthwhile 
to  consider  other,  if  related  schemes.  The  following  suggest  themselves: 

(1)  If  a  conventional  field  emitter  is  placed  roughly  one  tip  radius  above 
a  flat  conducting  surface,  the  field  on  the  plane  immediately  under  the  emit¬ 
ter  apex  is  quite  high  but  falls  off  rapid ly  with  the  distance  along  the  plane 
Thus  a  relatively  small  area  is  "illuminated"  and  conventional  field  emission 
may  be  obtained  from  it.  If  the  flat  surface  is  covered  with  adsorbate  which 
is  mobile,  a  "conventional"  field  emission  current  fluctuation  correlation 
function  can  be  obtained.  (2)  If  a  sufficiently  finely  etched  cylindrical 
wire  is  placed  near  and  parallel  to  the  surface  a  long  rectangular  region 
of  high  field  will  be  created  and  diffusion  across  the  narrow  dimension  can 
be  observed  as  a  one-dimensional  correlation  function  [3].  If  the  wire  is 
then  rotated  90°  the  other  component  of  the  diffusion  tensor  can  be  determined 
These  techniques  differ  from  our  present  procedures  only  in  that  a  flat 
substrate  is  the  field  emitter  and  a  highly  curved  object  serves  to  generate 
the  requisite  fields  over  relatively  small  areas.  There  is  no  insuperable 
problem,  given  the  state  of  the  STM  art,  in  implementing  the  first  scheme. 

It  may  be  possible  to  implement  the  second  by  etching  a  very  fine  loop  so 
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that  its  end  approximates  a  flat  cylinder,  or  by  etching  a  blade- like 
emitter  to  simulate  a  short,  thin  half-cylinder.  The  main  advantages  of 
these  schemes  are  that  the  states  of  the  tip  or  cylindrical  wire  are 
relatively  unimportant  since  they  only  serve  to  generate  a  field  and  to 
accept  electrons.  It  will  turn  out,  however,  that  the  requirements  on 
vibrational  isolation  and  drift  stability  are  still  very  severe  and  in  par¬ 
ticular  require  that  the  cylinder  be  short  in  length.  We  now  discuss  the 
forms  of  the  current  correlation  functions  f^  to  be  expected  for  these  two 
configurations . 


Tip  Above  a  Plane 


The  envisaged  geometry  is  shown  in  Fig.  1.  For  present  purposes  we 
assume  that  the  tip  can  be  replaced  by  a  conducting  sphere.  The  field  at 
the  plane  surface  is  then  given  by 

F(r)/Fo  =  (1  +  (r/zo)V3/2  (4) 

z  is  the  distance  of  the  sphere's  center  from  the  surface,  r  is  the 
o 

radial  distance  along  the  surface  from  the  origin  shown  in  Fig.  1;  Fq  is 
the  field  on  the  surface  at  r  =  0  and  is  given  by 


r  (2z  -  r  ) 

Fq  “  -v  - -  (5) 

z  (z  -  rj 
o  o  t 

where  V  is  the  applied  voltage  and  rt  the  tip  radius.  It  is  also  useful  to 
know  the  ratio  of  the  plane  surface  to  tip  apex  field  : 


/F  .  2  ♦  1 


tip  o 


( 2a- 1) 


where  a  =  zQ/rt.  If  a  =  2  this  ratio  is  2  so  that  the  apex  field  will  not  be 
high  enough  to  cause  field  desorption,  field  evaporation  or  field  ionization  if 

7 

the  field  at  the  plane  surface  is  3-6  x  10  v/cm,  i.e.  typical  for  field  emission 


V  -  •'*  A  A  .*•  .*• 


We  consider  next  the  correlation  function.  Since  relatively  large 


areas  of  emission  will  be  involved  we  work  out  only  the  analogue  of  gj(t) 
in  Refs.  1  and  2.  We  then  have  from  the  Fowler-Nordheim  equation 


In  j  =  InB  -  6.8  x  10 7  4>3/2/F 
6  In  j  =  Cj  6n  +  c^  6<t> 

=  (Cj  +  2tt  Pc^dn 


where  n  is  adsorbate  density  in  atoms/cm  , 

3lnB 
C1  ~  9n 

c2  =  |  6.8  x  10 7  «J»^/F 
and  2ttP  = 


C7) 

(8) 


(9) 

(10) 
(ID 


Here  4>  is  work  function,  P  adsorbate  dipole  moment,  B  a  preexponential,  and  j 
current  density.  Since  <Slnj  ^  6j/j  we  have  from  Eqs  .  8  and  4 


<Ai(o)Ai(t)>  =  /~d2r/“d2r' 
o  o 

=  (c,  +  2TTPC  )2  <j  >2  /“ d2 

1  l  O  o 


j(r)j(r")  <6j(r,o)  j(r  ,t) 

-3/2  a[(r/zo)2  ♦  (r'/t  ) 

r  /  d  r  e 


2] 

<6n(r,o)  n(r',t)> 


(12) 


Here  i  is  the  emitted  total  current,  Ai  its  fluctuation,  j  =  <j>  mean  current 

*  Jo  Jo 

density  at  r  =  0  and 


7  3/2 

a  =  6.8  x  10  «J»  /F  (13) 

We  have  used  the  fact  that  j  decreases  very  rapidly  with  decreasing  F  to 
expand  F/Fq  in  Eq .  4.  Since  the  term  exp  -a(r/zQ)2  still  falls  off  very 
rapidly  with  increasing  r  retaining  the  upper  limit  as  00  in  Eq.  12  introduces 
only  a  very  small  error.  The  quantity  <6n(r,o)6n(r',t)>  is 


S  e"lr"r"|2/4Dt 

o _ 


<6(r,o)6n(r'',t)>  = 


(14) 


(15) 


S  =  r  =  k  T  K  <n> 

o  A  B 


an  intensive  quantity  [2,7].  Here  <(6N)  >  is  the  total  mean  square  number 
(not  density)  fluctuation  in  an  area  A  and  K  the  adlayer  compressibility. 
The  total  current  <i>  is  found  from 


<i>  =  2it  /  rdrj  e 
o  ° 


-ot(l+(r/zo)  2)  3/2 


(16) 


2  2  - 


•=—  Z  TTJ 
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so  that  the  current  correlation  function  f. 

l 


f.(t)  ;  <Ai(o)M(t)> 


C17) 


<1> 


is  found  from  Eq.  11,  after  integration  over  the  angle  variables  [l]  to  be 

2 


sQ(ci  +  ^irPc-)  oo  oo  _ ( i+t” r 

ft(t)  =  - -  f  pdp  l  p'dpI^(2pp7T)  e  )(P  P  ) 


ff  Z  T 
O 


O  O 


Here 


■(£)■ 


(19) 


p  =  r/z' 


T  =  t/T 


T  =  z"*74D 
O  O 


(20) 

(21) 

(22) 


Iq  is  the  Bessel  function  Jq  of  imaginary  argument.  By  further  letting 


p/x  =  x  we  have  the  integral 


T  A dx  A'dx'I  C2xx')d-(1  *  T)fx2  *  x'2) 
o  o  o 


(23) 


We  multiply  the  integrand  by  JQ  (yx')  and  let  y  0  after  integration  over  x 


which  then  yields,  since  Jq(0)  =  1, 


lim  /  x'dx'  e~ex  I  (2xx')  J  (yx~) 


(24) 


4X-V 

, .  1  4e 

=  iim  ^  e  JQ 


ft1) 


(24) 


1  1  +  T 

e 


2(1  ♦  x) 

where  e  =  1  +  x.  Thus  the  integral  over  x  becomes  trivial  and 
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S  (c  +  2tiPc  ) 

f(t/xj  =  - ^ - 


1 


X  o' 


TT  Z 


8  ♦  4(t/xo) 


(25) 


A  comparison  of  f^(t/xQ)  according  to  Eq.  25  with  gj(t/xQ)  for  a  circular 
area  of  radius  rQ  is  shown  in  F'g.  2.  Clearly  the  functions  are  quite  similar 
but  not  identical.  It  is  interesting  that  the  long  time  behavior  of  f^  is 
the  expected  one  [2]  namely 


fi(t)  - 


Sq(Ci  +  2ttPc2)  4D 

IT  t 


(26) 


but  that  it  does  not  contain  zq,  so  that  in  principle,  D  can  be  determined  with¬ 
out  knowing  z  . 


Cylinder  Above  A  Plane 

It  is  assumed  that  the  radius  of  the  cylinder  can  be  rQ  'v  2-3  x  10  5  cm 

so  that  its  macroscopic  length,  which  we  will  call  2b  can  be  considered  infinite, 

for  purposes  of  estimating  the  field  on  the  surface  underneath  the  cylinder, 
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even  if  b  =  10  -  10  cm.  With  this  assumption  the  field  on  the  surface  at 

x  =  0  (notation  is  shown  in  Fig.  3)  is 

F  =  -2V/[z  ln(2(z  /r  )-l)]  (27) 

O  O  O  U 

and  the  ratio 

F(x)/Fq  =  (1+  (x/zo)2)-1  (28) 

The  ratio  of  the  field  at  the  cylinder  surface  (for  x  =  0)  to  that  on  the 
plane  is 


F  ,/F 


2 


(29) 


where  a  =  r  /z  ,  r  being  the  cylinder  radius.  If  a  =0.5  the  field  at 

O  O  0 

the  cylinder  surface  is  2.7  Fq  so  that  again  field  desorption  etc.,  will 
not  take  place.  The  total  current  field  emitted  from  the  substrate  is 

i  =  jQ  2b  (tt/cx)^  zo  (30) 


Here  j  is  the  current  density  corresponding  to  Fq  and  a  is  defined  by  Eq.  13. 

If  the  cylinder  is  parallel  to  one  of  the  principal  axes  of  the  diffusion 
tensor  the  density  (and  current)  correlation  function  decomposes  into  a  pro¬ 
duct  of  two  1-dimensional  functions.  We  then  have,  by  analogy  to  the  previous 
development  in  this  paper 


„  So(Cl  *  2lrPC2>2  2  *  (*'/S)2ye-|x-X'1  /4D**t 

1  =  4b  V/CO  x/“ x  e 

O  *  V  \ 


/b 


'|y-y"l2/4Dyyt 

4tFd  t 

yy 


The  integrals  over  y ,  y'  yield  2b  since  in  the  time  of  an  experiment  fy(t)  'v- 
f  (0)  [3],  The  quadratures  over  x,  x'  are  elementary  and  one  obtains  finally 


S  (c,  +  2ttPc_) 

f  ftl  -  u  1  L  *  f 

i  V  7  V 

2b(n/a)  2  zq 

it  x  +  2 

where 

T  =  t/x 

o 

(33) 

and 

2 

o  #  1 

To  =  a  D 

XX 

(34) 

The 

form  of  f^  is  compared  with 

that  for  a  rectangle  in  Fig.  4.  It  is 

apparent 

that  the  curves  are  remarkably  similar. 

The  long  time  behavior  of 

f^  is  given  by 

S  (c  +  2ttPc  )2 
MO  =  - r - — 

A  ;^» 

(l/2u)  (35) 
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which  again,  does  not  contain  z  .  It  is  interesting  to  note  that  f.  is 

o  x 

inversely  proportional  to  2b,  the  length  of  the  cylinder.  However,  the 
actual  current  fluctuations,  i.e.  <Ai(o)Ai(t)>  increase  with  b. 

<Ai(o)Ai(t)>  =  +  2ttPC2)  2< j  >*"  C 27T/ot)  2  zQb  (36) 

Thus  an  a.c.  measurement,  which  is  what  is  normally  performed  will  actually 
yield  a  bigger  signal  the  longer  the  cylinder,  although  stability  considera¬ 
tions  will  show  that  b  should  be  kept  as  small  as  possible. 

Determination  of  a  and  z 
_ o 

The  determination  of  activation  energies  requires  only  that  zq,  the  distance 
of  the  tip  or  cylindrical  wire  center  from  the  plane  substate  be  kept  constant 
as  the  substrate  temperature  is  varied.  A  determination  of  prefactors,  i .e. 
absolute  values  of  D  requires,  however,  that  a  and  zq  be  known  or  at  least 
approximated.  We  now  indicate  how  this  could  be  done: 

(1)  Tip  above  Plane 

For  this  case  the  current  emitted  by  the  plane  substrate  is  given  by 
Eq.  16,  which  can  be  written,  in  view  of  the  definition  of  a  by  Eq.  13,  as 


CkV 

6.8  x  107  4>3/2 


-6.8xl07c}>3/2/kV 


where 


k  =  F  /V 
o 


and  C  is  the  (image  corrected)  Fow ler-Nordheim  preexponential  for  a  clean 
metal  surface  [8] 

C  =  6.2*106^V-40','3/2  (39) 

p  +  4)  > 

The  exponential  term  in  the  RHS  of  Eq.  39  comes  about  because  the  image 


correction  term  v(y)  where 


y  =  3. 8xl0*4  F^/4> 
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which  multiplies  the  Fowler-Nordheim  exponent  turns  out  to  have  the  form 

v(y)  =  1  -  0.059  x  10'7  Fq  (41) 

so  that  the  image  correction  appears  as  a  factor  in  the  Fowler-Nordheim 
p re exponential . 

3 

It  is  therefore  possible,  at  fixed  zq  to  determine  lni/V  vs.  1/V 
(which  is  the  analogue  of  a  Fowler-Nordheim  plot  used  in  ordinary  field 
emission)  and  from  its  slope  to  determine  k,  assuming  that  0  is  known, 
for  the  clean  surface.  Thus  a  is  found  from  Eqs.  13  and  38  for  given  V. 
zq  can  then  be  found  from  the  total  current  and  Eqs.  37  and  39. 

It  is  also  possible  to  utilize  Eq.  5  since 

k  =  - 2 — L_  (42) 

z  (z  -r  ) 
ov  o  t 

if  the  emitter  radius  r^  is  known,  either  from  field  emission  measurements  from 
the  tip  at  large  distance  from  the  substrate  or  from  independent  determination. 
However,  Eq.  5  is  based  on  an  assumed  sphere-plane  geometry  and  probably  less 
accurate  that  the  method  described  previously  which  uses  only  emission  from 
the  substrate. 

(2)  Cylinder  Above  Plane 

For  this  case  a  and  k  =  Fq/V  can  be  determined  analogously  to  the  tip-plane 
5/2 

case  by  plotting  lni/V  vs.  1/V  as  can  be  seen  from  Eq.  30.  A  procedure 
analogous  to  that  outlined  for  the  tip-plane  case  could  then  be  used  but  turns 
out  to  be  largely  unnecessary  since 


k  . -L 

z 

o 


In  2(z  /r  )-l 
o  o' 


since,  for  appreciable  emission  the  logarithmic  terra  is  close  to  unity. 


•  A  w'*  ,**  .*•  ',*•  i  ••  »' *>*%•  «*•  *'*»  /*  ,  ^  »*•  »'*  .l’  «'»  ,**  4%  .**  %*•  „*«•*%  ,*• . ^  »**  AV.*  '•  . 
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Stability  Requirements 

It  is  necessary  also  to  examine  the  required  stability  of  tip  or  cylinder  both 
with  respect  to  vibrations  and  "long-time"  drift  since  small  changes  in  tip-substrate 
or  cylinder-substrate  distance  will  lead  to  changes  in  fields  and  hence  cur¬ 
rents,  which  could  interfere  with  the  diffusion  induced  fluctuation  signal. 

(1)  Tip  Above  Plane 
(a)  Vibrations 

We  ignore  the  effect  of  field  on  the  preexponential  term  of  the  cur¬ 
rent  and  write,  using  Eqs .  7  and  16 

■Slni  =  y  =  6.8  x  107  4>3/2  6(1/Fq) 


Use  of  Eq.  5  then  yields 


6i  _  *  2  3  a2-l  | 

T  6  2TT  ‘  “1 - 

L  *  -a  J 


where  a  =  z  /r  ,  so  that 
o  t 


-o  ■  m)  *  - 


We  must  now  compare  this  signal  with  /f^('o)  to  get  a  sense  of  its  importance. 
From  Eq.  25  we  have  that 


ft(o)  - 


C<5i)  2>  _  <(6N)  2>  (cl  +  27TPc2) 


c  2  2 
8  TT  -=—  Z 

3a  o 


We  now  assume  roughly  [7]  that 


<(6N)  >  -v  0.1  <N> 


=  0.1  A  <n> 


so  that 


f  2V 

I  0.1  <n>  (c.  +  2ttPcJ  1 
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where  we  have  used  a  =  z  /r  . 

o  t 

We  now  require  that 


Pf)  =  y  /T(o)  (50) 

\  /vibration 

where  y  is  a  tolerance  factor  whose  value  depends  on  the  details  of  the  vibra¬ 
tional  noise  and  the  dynamic  range  of  the  correlator.  If  the  noise  is  random 
and  the  dynamic  range  of  the  correlator  sufficient  y  >  10-50  might  still  allow 
extraction  of  a  diffusion  correlation  function  if  the  data  collection  time  is 
increased  accordingly.  If  vibrations  are  coherent  y  might  have  to  be  as  low 
as  0.1.  Combining  Eqs .  50,  49  and  46  we  find  assuming  a  =  2 


6z  -  y  — — 

°  ifiT 


+  2irPc2) 


15  2 

To  get  some  feeling  for  numbers  let  us  assume  <n>  n.  10  atoms/cm  and  a  =  20 

~  - 15  2 

as  before.  We  choose  ^  =  dlnB/dn  =  -2x10  cm  /atom  and  2itP  =  d<f>/dn  = 

lxio'15  eV  cm2/atom;  c2  =  10  (p^/F  if  F  is  expressed  in  volts/A;  we  let 
-1 

c 2  =  7.5  (eV)  corresponding  to  4>  =  5  eV  and  F  =  0.3  volt/A.  These  values 
correspond  approximately  to  oxygen  on  a  W(110)  plane  [9].  We  then  find  that 

o 

6zq  =  y  0 . 2  A.  If  y  ^  1  is  tolerable  the  vibrational  amplitude  is  within 
the  limits  of  STM  technology.  It  is  interesting,  incidentally,  that  6zq  is 
independent  of  zq  (or  r  )  and  depends  only  on  the  ratio  a,  although  not  dra¬ 
matically.  This  relative  independence  of  6zq  on  collector  radius  suggests  a 
way  of  circumventing  vibrational  noise  which  generally  is  of  low  frequency. 

If  the  collector  is  made  fairly  sharp  and  or  if  the  measurements  are  carried 
out  at  high  enough  temperature  to  make  tq  so  small  that  the  relevant  frequency 
domain  lies  outside  that  of  the  mechanical  vibrations  y  could  be  much  larger 
than  suggested  above. 

(b)  Drift 


In  addition  to  vibrations  one  might  also  be  faced  with  a  more  or  less 
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steady  drift.  We  will  assume  that  Btq  is  the  minimum  collection  time,  i.e. 

that  B  runs  of  duration  tq  each  must  be  taken  to  obtain  a  meaningful  average. 
3 

B  is  probably  ^  10  depending  on  random  noise.  We  will  also  assume  that 
in  time  B  tq  the  mean  value  of  the  emission  current  must  change  by  <0.1  <i>. 
We  then  have,  combining  Eqs.  45  and  22  and  assuming  a  =  2,  that 
dz 


o  <0.13  D 


dt  it  B  r^ 


(52) 


If  r  =  300  A,  D  =  lO-11  cm2sec_1  and  B  =  103  we  find  dz  /dt  =  0.01  A/sec. 
t  o 

While  this  seems  very  stringent  it  should  be  remembered  that  it  is  possible  to 
regulate  zq  at  rates  (i.e.  frequencies)  which  do  not  interfere  with  the  correla¬ 
tion  signal  by  adjusting  the  servomechanism  to  hold  the  mean  current  fixed, 
even  though  the  a.c.  component  is  being  used  to  obtain  a  correlation  signal. 


(2)  Cylinder  Above  Plane 
(a)  Vibrations 


For  this  case  one  finds  from  Eqs.  7  and  27  that 

Si/i  =  a  dz  /z_ 
o  o 


(53) 


since  the  change  in  the  logarithmic  term  in  Eq.  27  can  be  ignored.  An  analysis 
similar  to  that  leading  to  Eq.  51  then  gives 


5zq  =  4.4x  10"9  Y(zo/b)’S  (54) 

assuming  the  same  numerical  values  of  a,  <n>,  c^,  and  2ttP  used 

previously.  If  y  =  10  we  find  that  6zq  =  1.42  A  if  ZQ/b  =0.1  and  0.45  A  if 

•  -3 

z/b  =  0.01:ifz  =  4000  A  b  could  therefore  be  as  much  as  4x10  cm  and 

o  o  — 

still  give  a  tolerable  value  of  6zq.  However,  for  cylindrical  geometry  tq  is 
necessarily  larger  than  for  the  hemispherical  collector,  since  we  have 
assumed  that  the  cylinder  radius  is  larger.  Thus  the  domain  of  frequencies 
of  relevance  to  the  correlation  function  will  be  lower  and  may  include  vibra¬ 
tional  frequencies,  so  that  y  may  have  to  be  <  10,  the  value  assumed  here. 


.Oil  — 


Nevertheless,  there  should  be  some  ranges  of  D  which  can  be  investigated 
with  the  cylindrical  collector,  at  least  as  far  as  vibration  amplitudes 
are  concerned. 


(b)  Drift 


We  use  Eq .  34  for  tq  and  find 

0.1  D 

dz  /dt  <  - - ** 

O  6  IT  2 


•  3  7  • 

If  we  assume  zq  =  4000  A  and  $  =  10  we  find  dzQ/dt  <  8x10  D  A/sec.  This  very 

low  value  puts  a  rather  definite  limit  on  measurable  D  values,  i.e.  >  10"9 
2  -1 

cm  sec  ,  hut  this  is  consistent  with  the  fact  that  tq  is  larger  than  for 
the  hemispherical  collector. 


Field  Effects 

We  consider  finally  the  effect  of  high  field,  not  on  diffusion  per  se 
but  on  the  changed  adsorbate  concentration  arising  from  field-dipole  and 
polarization  effects.  In  a  region  of  high  field  the  adsorbate  density  is 

n(f)  =  n(0)  exp(Js  a  F2/kT  +  P.F/kT)  (56) 


where  n(F)  and  n(0)  are  the  adsorbate  densities  in  the  field  and  field  free 
regions  respectively;  a  is  the  polarizability  and  P  the  dipole  moment.  If 
the  field  is  appropriate  for  field  emission  and  the  negative  end  of  the 
dipole  moment  points  outward,  i.e.  if  the  work  function  is  raised  by  adsorption 
P.F  in  Eq.  56  is  positive;  if  there  is  a  work  function  decrease  caused  by 
adsorption  P.F  in  Eq .  56  is  negative.  For  P  in  Debye,  a  in  A3  and  F  in  volts/ 

A  Eq.  56  takes  the  form 

n(f)  =  n(0)  exp[(403  aF2  +  2.42  x  103  P.F)/T]  (57) 

.5  is  7 

For  a=1A,P  =  0.27  Debye  (corresponding  to  Ad)  =  1  eV  for  n  =  10  atoms /cm  ) 

and  F  =  0.4  volts/A  we  find,  at  300  K  n(f)/n(0)  =  2.96  and  1.7  at  600  K  This 


is  a  substantial 
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increase  and  indicates  that  the  concentration  in  the  probed  region  may  be 
considerably  enhanced.  This  is  no  fundamental  obstacle  since  the  concen¬ 
tration  in  the  probed  region  can  be  measured  by  determining  the  work  func¬ 
tion  change. 

Diffusion  itself  will  also  be  affected  in  two  ways.  First  D  is 
generally  a  function  of  n  and  thus  having  different  values  of  n  in  the 
probed  region  and  beyond  may  affect  measurements.  Second,  roughly  speaking 
atoms  now  undergo  biassed  random  walks  since  there  is  a  field  contribution 
on  to  their  binding  energy.  Nevertheless,  the  situation  is  not  as  bad  as 
it  might  appear,  for  roughly  the  same  reason  why  field  variation  is  not  a 
serious  problem  in  the  fluctuation  method  employing  a  field  emitter  as  sub¬ 
strate:  In  that  case  the  field  also  drops  off  as  one  moves  away  from  the 
tip  apex  but  the  region  over  which  diffusion  is  of  importance  is  confined  to 
that  region  of  the  tip  where  the  field  is  nearly  uniform.  In  the  present 
case  field  emission  is  a  much  more  sensitive  function  of  F  and  hence  r  than 
the  contribution  of  F  to  binding  energy;  in  particular  we  have  picked  a 

rather  high  polarizability  in  the  example  and  in  most  cases  of  interest 
2 

!jaF  can  probably  be  neglected,  so  that  only  the  P.F  term  remains.  Then 
field  emission  occurs  from  a  relatively  small  region  of  almost  uniform  F 
while  diffusion  will  occur  with  nearly  constant  D  from  a  much  larger  region 

with  nearly  constant  n.  It  is  easy  to  show,  for  the  tip  above  plane 

case,  from  Eqs.  4  and  57  for  the  values  of  a,  P  and  F  used  that  n/n(0) 

2 

will  have  decreased  from  2.96  at  (r/z  )  =  0  to  2.27,  i.e.  by  only  30% 

o 

2 

at  (r/z^  )  =5,  which  corresponds  to  2.24  times  the  radius  of  the  probed 

region.  Nevertheless,  the  field  effect  cannot  be  neglected  and  must  be 
considered  case  by  case.  For  the  single  site  method  where  very  small 
fields  can  in  principle  be  used  it  can  probably  be  neglected  altogether. 

It  should  perhaps  be  pointed  out  that  the  effect  we  arc  considering  here 


is  different  from  that  of  a  uniform  high  field  everywhere  on  the  activation 

energy  of  diffusion.  This  latter  effect  concerns  the  differences  in  P,  F 
2 

and  %  F  at  the  potential  minima  and  saddle  points  of  diffusing  atoms  and 
is  important  only  for  high  dipole  moments  and  very  high  polarizabilities, 
such  as  are  encountered  for  adsorbed  alkali  atoms. 

Cone lusion 

The  foregoing  has  indicated  that  it  may  be  possible  to  measure  surface 
diffusion  coefficients  of  adsorbates  by  extensions  of  the  field  emission 
current  fluctuation  correlation  method.  Three  possible  schemes  have  been 
proposed  and  analyzed.  Each  is  most  suitable  for  a  different  range  of  dif¬ 
fusion  coefficients,  since  the  relaxation  times  involved  are  different.  Each 
method  has  somewhat  different  constraints  imposed  by  vibrational  noise  and 
drift.  In  general  vibrational  noise  and  drift  within  a  given  scheme  are 
antagonistic:  In  order  to  overcome  vibrational  noise  longer  data  collection 
times  will  be  required  and  this  requires  less  drift  or  better  drift  compen¬ 
sation.  However,  it  appears  that  each  method  could  actually  be  used  over  a 
limited  range  of  D  values. 

No  attention  has  been  paid  to  actual  design  considerations,  i.e.  the 
need  to  prepare  clean  substrates,  maintain  ultrahigh  vacuum  and  control 
coverage,  for  instance  by  work  function  measurements, or  how  to  achieve  tempera¬ 
ture  variation  without  appreciably  changing  the  collector-substrate  distance. 
These  problems,  although  not  trivial,  seem  soluble.  For  instance,  emission 
from  the  collector  can  always  be  used  to  control  its  distance  from  the 
substrate,  even  when  the  latter  is  a  semiconductor  with  temperature  dependent 
field  emission.  Thus  there  is  some  hope  that  the  schemes  proposed  here  can 
actually  be  implemented. 
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Figure  Captions 

1)  Schematic  diagram  of  tip-plane  geometry:  r  radius  of  (assumed  spherical) 
tip;  r  vector  in  the  surface  from  origin  directly  under  the  tip  apex  in 
the  plane;  zq  distance  from  emitter  center  to  surface. 

2)  Comparison  of  correlation  function  gjCt/x  )  for  a  circular  aperature  of 

2 

radius  r  ,  for  which  t  =  r  /4D  (solid  line  and  scale  marked  (a))  with 
o  o  o 

normalized  correlation  function  f^  =  2/(2  +  t/xQ)  for  the  tip-above 

plane  geometry  (points  and  scale  marked  (b)) .  For  this  latter  case 
2 

Xq  =  (2/3a)zQ/4D  and  the  horizontal  axis  has  been  displaced  for  best  fit 
with  the  g1(t/xQ)  curve. 

3)  Schematic  diagram  of  the  cylinder-plane  geometry:  rQ,  cylinder  radius; 
zq,  distance  from  cylinder  axis  to  surface;  x-y,  coordinate  frame  in  the 
surface  plane;  2b,  length  of  cylinder. 

! 

4)  Comparison  of  1-dimensional  correlation  function  for  the  side  of  a 

rectangle  of  width  2a  (solid  line  and  scale  marked  (a))  for  which 

2  i< 

x  =  a  /D  with  normalized  correlation  function  [2/(2  +  4irt/x  )]J 
o  xx  o 

2 

(points  and  scale  marked  (b) ) ,  for  which  xq  =  nzo/aDxx.  The  horizontal 
axis  has  been  displaced  for  best  fit  of  points  to  line. 
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